function new=Gauss1DFilter(image,Sigma)
maskSize=2 * int32(3.7 * Sigma - 0.5) + 1;
GF1D = fspecial('gaussian', [1 double(maskSize)], Sigma) ;
horBlur = imfilter(image, GF1D) ;
new = imfilter(horBlur, GF1D') ;